Excess noise for coherent radiation propagating through amplifying random media 
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A general theory is presented for the photodetection statistics of coherent radiation that has been 
amplified by a disordered medium. The beating of the coherent radiation with the spontaneous 
emission increases the noise above the shot-noise level. The excess noise is expressed in terms of the 
transmission and reflection matrices of the medium, and evaluated using the methods of random- 
matrix theory. Inter-mode scattering between N propagating modes increases the noise figure by 
up to a factor of N, as one approaches the laser threshold. Results are contrasted with those for an 
absorbing medium. 

PACS numbers: 42.50.Ar, 42.25.Bs, 42. 25. Kb, 42.50.Lc 



I. INTRODUCTION 

The coherent radiation emitted by a laser has a noise 
spectral density P equal to the time-averaged photocur- 
rent I. This noise is called photon shot noise, by anal- 
ogy with electronic shot noise in vacuum tubes. If the 
radiation is passed through an amplifying medium, P in- 
creases more than I because of the excess noise due to 
spontaneous emission |Q. For an ideal linear amplifier, 
the (squared) signal-to- noise ratio I 2 /P drops by a fac- 
tor of two as one increases the gain. One says that the 
amplifier has a noise figure of 2. This is a lower bound 
on the excess noise for a linear amplifier 0. 

Most calculations of the excess noise assume that 
the amplification occurs in a single propagating mode. 
(Recent examples include work by Loudon and his 
group . ) The minimal noise figure of 2 refers to this 
case. Generalisation to amplification in a multi-mode 
waveguide is straightforward if there is no scattering be- 
tween the modes. The recent interest in amplifying ran- 
dom media || calls for an extension of the theory of 
excess noise to include inter-mode scattering. Here we 
present such an extension. 

Our central result is an expression for the probability 
distribution of the photocount in terms of the transmis- 
sion and reflection matrices t and r of the multi-mode 
waveguide. (The noise power P is determined by the 
variance of this distribution.) Single-mode results in the 
literature are recovered for scalar t and r. In the ab- 
sence of any incident radiation our expression reduces to 
the known photocount distribution for amplified sponta- 
neous emission ||. We find that inter- mode scattering 
strongly increases the excess noise, resulting in a noise 
figure that is much larger than 2. 

We present explicit calculations for two types of ge- 
ometries, waveguide and cavity, distinguishing between 
photodetection in transmission and in reflection. We 
also discuss the parallel with absorbing media. We use 
the method of random-matrix theory to obtain the 
required information on the statistical properties of the 
transmission and reflection matrices of an ensemble of 



random media. Simple analytical results follow if the 
number of modes N is large (i.e. for high-dimensional 
matrices). Close to the laser threshold, the noise fig- 
ure T exhibits large sample-to-sample fluctuations, such 
that the ensemble average diverges. We compute for ar- 
bitrary N > 2 the distribution p(T) of T in the ensemble 
of disordered cavities, and show that T — N is the most 
probable value. This is the generalisation to multi-mode 
random media of the single-mode result T = 2 in the 
literature. 



II. FORMULATION OF THE PROBLEM 

We consider an amplifying disordered medium em- 
bedded in a waveguide that supports N(u>) propagat- 
ing modes at frequency u> (see Fig. |l|). The amplifica- 
tion could be due to stimulated emission by an inverted 
atomic population or to stimulated Raman scattering [jy . 
A negative temperature T < describes the degree of 
population inversion in the first case or the density of 
the material excitation in the second case || . A complete 
population inversion or vanishing density corresponds to 
the limit T — *■ from below. The minimal noise figure 
mentioned in the introduction is reached in this limit. 
The amplification rate l/r n is obtained from the (neg- 
ative) imaginary part e" of the (relative) dielectric con- 
stant, 1/T a = w|e"|. Disorder causes multiple scattering 
with rate l/r s and (transport) mean free path I = ct s 
(with c the velocity of light in the medium) . We assume 
that t s and r a are both 3> 1/w, so that scattering as 
well as amplification occur on length scales large com- 
pared to the wavelength. The waveguide is illuminated 




FIG. 1. Coherent light (thick arrow) is incident on an am- 
plifying medium (shaded), embedded in a waveguide. The 
transmitted radiation is measured by a photodetector. 
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from one end by monochromatic radiation (frequency too, 
mean photocurrent Iq) in a coherent state. For simplicity, 
we assume that the illumination is in a single propagating 
mode (labelled mo). At the other end of the waveguide, 
a photodetector detects the outcoming radiation. We as- 
sume, again for simplicity, that all N outgoing modes are 
detected with equal efficiency a. 

We denote by p(n) the probability to count n photons 
within a time r. Its first two moments determine the 
mean photocurrent / and the noise power P, according 
to 
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P 



lim — ^ 

r — >oo T 
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(2.1) 



(The definition 

p = j^dtsimm 



of P is equivalent to 
with SI = I — I the fluctuating 



part of the photocurrent.) It is convenient to compute 
the generating function F(£) for the factorial cumulants 
Kj, defined by 



F(0 



(2.2) 



One has n — Ki, n 2 = K2 + «i(l + Ki). 

The outgoing radiation in mode n is described by an 
annihilation operator a° ut (w), using the convention that 
modes 1,2, ... ,N are on the left-hand-side of the medium 
and modes iV+1, . . . , 2N are on the right-hand-side. The 
vector a° ut consists of the operators a° ut , og"*, . . . , a^jy. 
Similarly, we define a vector a m for incoming radiation. 
These two sets of operators each satisfy the bosonic com- 
mutation relations 

[a n (u>), <4i,(w')] = 6 nm 5(uj - u>') , (2.3a) 
[a n (cj),a m (cj')] = 0, (2.3b) 

and are related by the input-output relations [HHH] 

a out (uj) = S(io)a ia (uj) + V(w)c + (w) . (2.4) 

We have introduced the 2N x 2N scattering matrix S, 
the 2N x 27V matrix V , and the vector c of 2N bosonic 
operators. The scattering matrix S can be decomposed 
into four TV x N reflection and transmission matrices, 



r' if 
t r 



(2.5) 



Reciprocity imposes the conditions t' = t T , r = r T , and 

/ IT 

r = r . 

The operators c account for spontaneous emission in 
the amplifying medium. They satisfy the bosonic com- 
mutations relation (2.3), which implies that 



VV^ = SS f - 1 
Their expectation values are 



(2.6) 



(c n M4>')) = -5 nm S(iu - u/)f(u, T) , (2.7) 
with the Bose-Einstein function 

f{uj,T) = [exp{hw/kT) - If 1 (2.8) 
evaluated at negative temperature T (< 0). 



III. CALCULATION OF THE GENERATING 
FUNCTION 

The probability p(n) t hat n photons are counted in a 
time r is given by M.O] 



p(n) = -(: W n e- W :) 



(3.1) 



where the colons denote normal ordering with respect to 
a ou \ and 



T 2N 

W = a dt £ a° ut+ (*K ut (i) 

*' n=N+l 



(t) = (27T)- 1 / 2 / dwe- Mt < ut H 



The generating function (2.2) becomes 
F(0 = ln(: e* w :) . 



(3.2) 
(3.3) 

(3.4) 



Expectation values of a normally ordered expression 
are readily computed using the optical equivalence theo- 
rem [jl2| . Application of this theorem to our problem con- 
sists in discretising the frequency in infinitesimally small 
steps of A (so that lu p — pA) and then replacing the an- 
nihilation operators a}ff(u> p ), c n {uj p ) by complex numbers 
a^ p , c np (or their complex conjugates for the correspond- 
ing creation operators). The coherent state of the inci- 
dent radiation corresponds to a non-fluctuating value of 
alf p with |ajf p | 2 = S nmo S PPo 2nI /A (with u = p A). The 
thermal state of the spontaneous emission corresponds 
to uncorrelated Gaussian distributions of the real and 
imaginary parts of the numbers c np , with zero mean and 
variance ((Rec„ p ) 2 ) = ((Imc np ) 2 ) = -\f(u p ,T). (Note 
that / < for T < 0.) To evaluate the characteris- 
tic function (3.4) we need to perform Gaussian averages. 
The calculation is described in the appendix. 

The result takes a simple form in the long-time regime 
lj c t 3> 1, where lo c is the frequency within which S(u>) 
does not vary appreciably. We find 



/ ln\\l - a£f(l - rri -tP)\\dw 

2tt Jo 
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F exc (0 = a^rlo 

x (V [1 - a£/(l - rr f 
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(3.5) 
(3.6) 
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denotes the determinant and (- - /momo 



where || • 

mo, mo-element of a matrix. In Eq. ( 
/, t, and r are to be evaluated at u 



jifiomo the 
i the functions 
wo- The inte- 
gral in Eq. (3.5) is the generating function for the pho- 
tocount due to amplified spontaneous emission obtained 
in Ref. ||. It is independent of the incident radiation 
and can be eliminated in a measurement by filtering the 
output through a narrow frequency window around uiq. 
The function F exc (£) describes the excess noise due to the 
beating of the coherent radiation with the spontaneous 
emission. The expression ( |3.6| ) is the central result of this 
paper. 

By expanding F(£) in powers off we obtain the fac- 
torial cumulants, in view of Eq. ( [2.2| ). We will in what 
follows consider only the contribution from -F'exc(f), as- 
suming that the contribution from the integral over u has 
been filtered out as mentioned above. We find 

n k = W^t/^/o^ (1 - rrt - ttt)*- 1 t ] momo , (3.7) 

where again u> = u>o is implied. The mean photocurrent 
I = Ki/r and the noise power P = {k% + Ki)/t become 



I = al (tH) 

P oxc = 2a 2 fI [S (1 



P = I + P c: 



tt f ) t], 



(3.8) 



The noise power P exceeds the shot noise I by the 
amount P exc - 

The formulas above are easily adapted to a measure- 
ment in reflection by making the exchange r — ► f', 
t — > r' . For example, the mean reflected photocurrent 
is I = alo (r''r') , while the excess noise is 



Pexc - 2a 2 fl [r'^ (a - r'r't - t't't) r% 



(3.9) 



IV. NOISE FIGURE 

The noise figure T is defined as the (squared) signal- 
to-noise ratio at the input Iq/Pq, divided by the signal- 
to-noise ratio at the output, I 2 /P. Since Po = Iq for 
coherent radiation at the input, one has J- = (P exc + 
T)Io/I 2 , hence 



T=-2f 



(tfrrH + tHtH) i 
(ttt) 



2af 



a (ttt) 



(4.1) 



The noise figure is independent of Iq. For large amplifi- 
cation the second term on the right-hand-side can be ne- 
glected relative to the first, and the noise figure becomes 
also independent of the detection efficiency a. The min- 
imal noise figure for given r and t is reached for an ideal 
detector (a = 1) and at complete population inversion 

(/ = -!)• 

Since (tWt + tttttt) momo = £ fe |(tt r ) moi 



E fe l(* f *) 



mo&i I 



> (* + *)m m > one has T ^ ~ 2 -f for lar S e 



amplificatio n (w hen the second term on the right-hand- 
side of Eq. (hi) can be neglected). The minimal noise 
figure T = 2 at complete population inversion is reached 
in the absence of reflection [(t^r) mg k — 0] and in the ab- 
sence of inter-mode scattering [(ttt) mo fc = if k ^ mg 
This is realised in the single-mode theories of Refs. [ 
Our result (4.1) generalises these theories to include scat- 
tering between the modes, as is relevant for a random 
medium. 

These formulas apply to detection in transmission. For 
detection in reflection one has instead 



T=-2f 



(r'H't'^r' + r't r'r't r >) l + 2af 



(r'Vr 



a (r'tr') r 



(4.2) 



Again, for large amplification the second term on 
the right-hand-side may be neglected relative to 
the first. The noise figure then becomes small- 
est in the absence of transmission, when T = 
-2/(r't r'r't r') momo (r't r') m 2 omo > _ 2/ . The minimal 
noise figure of 2 at complete population inversion requires 
(r't r'r't r') momo = (r'V)^ oTOo , which is possible only in 
the absence of inter-mode scattering. 

To make analytical progress in the evaluation of T we 
will consider an ensemble of random media, with different 
realisations of the disorder. For large N and away from 
the laser threshold, the sample-to-sam ple fluctu ation s in 
numerators and denominators of Eqs. ([0]) and (4.2) are 
small, so we may average them separately. Furthermore, 
the "equivalent channel approximation" is accurate for 
random media which says that the ensemble aver- 
ages are independent of the mode index mo. Summing 
over mo, we may therefore write T as the ratio of traces, 
so the noise figure for a measurement in transmission be- 
comes 



-2fN 



-< tr (ttrr f t + t^ttf) y 
-< trttf >- 2 



■N 



l + 2af 
a -< tvtU y ' 
(4.3) 



and similarly for a measurement in reflection. The brack- 
ets -<•••>- denote the ensemble average. 



V. APPLICATIONS 

A. Amplifying disordered waveguide 

As a first example, we consider a weakly amplifying, 
strongly disordered waveguide of length L (see inset of 
Fig. g) . Averages of the moments of rrt and ttt f or ^ n j s 
system have been computed by Brouwer [Q as a function 
of the number of propagating modes N, the mean free 
path I, and the amplification length f a = \fDr a , where 
l/r a is the amplification rate and D = d/3 is the diffu- 
sion constant. It is assumed that 1/N <C i/£ a "C 1 but 
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FIG. 2. Noise figure of an amplifying disordered waveguide 
(length L, amplification length £ Q ) measured in transmission 
(solid line) and in refl ecti o n (d ashed line). The curves are 
computed from Eqs. ( |}.l[ )-( [5^ ) for a = 1, f = —1, and 
L/l = 10. The laser threshold is at L/£ a — tt. 



the ratio L/£ a = s is arbitrary. In this regime, sample- 
to-sample fluctuations are small, so the ensemble average 
is representative of a single system. 

The results for a measurement in transmission are 



Aal s 

SL sins 

2a 2 l 
Pcxc = -7TT fhs 



3L 



+ 



sm s 
s cotan s — 1 



2s — cotan s 
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sm s 



sin 3 s 



• 4 

sm s 



For a measurement in reflection, one finds 
Al 



I = 
p — 

± cxc — 



al 

2a 2 l 
3L 



1 



3L 



-s cotan s 



2 cotan s 



1 



s cotan s — 1 
sin 3 s 



sm s 
s 



cotan s 
sm s 



• 4 

sm s 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



The noise figure J- follows from T = (P CX c + I)Io/I 2 - It 
is plotted in Fig. ||. One notices a strong increase in T 
on approaching the laser threshold at s = 7r. 



B. Amplifying disordered cavity 

Our second example is an optical cavity filled with an 
amplifying random medium (see inset of Fig. ^). The 
radiation leaves the cavity through a waveguide support- 
ing N modes. The formulas for a measurement in reflec- 
tion apply with t = Q because there is no transmission. 
The distribution of the eigenvalues of r^r is known in 
the large- N limit Eq| as a function of the dimensionless 



FIG. 3. Noise figure of an amplifying disordered cavity, 
connected to a photodetector via an iV-mode waveguide. The 
curve is the result (|5.9|), as a function of the dimensionless 
amplification rate 7. (Ideal detection efficiency, a = 1, and 
full population inversion, / = —1, are assumed in this plot.) 
The laser threshold occurs at 7 — 1. 



amplification rate 7 = 2-7T /Nr a \u> (with Aw the spacing 
of the cavity modes near frequency u>q). The first two 
moments of this distribution are 



iV _1 -<! trrV y = 



1 



1-7 

»r 1 t + 2 7 2 - 2 7 + 1 

(1 - 7 ) 4 

The resulting photocurrent has mean and variance 
1 



I = al 



1-7 



P cxc = 2a 2 // o7 



7 - 7 2 - 1 



(1-7) 4 ' 

The resulting noise figure for a — 1 and / = — 1, 
1 - 7 + 7 2 + 7 3 



(1-7) 



(5.5) 
(5.6) 

(5.7) 
(5.8) 

(5.9) 



is plotted in Fig. |[ Again, we see a strong increase of T 
on approaching the laser threshold at 7 = 1. 



VI. NEAR THE LASER THRESHOLD 

In the previous section we have taken the large- N limit . 
In that limit the noise figure diverges on approaching the 
laser threshold. In this section we consider the vicinity 
of the laser threshold for arbitrary N. 

The scattering matrix S(lo) has poles in the lower half 
of the complex plane. With increasing amplification, the 
poles shift upwards. The laser threshold is reached when 
a pole reaches the real axis, say at resonance frequency 
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LUth ■ For w near uj t h the scattering matrix has the generic 
form 



5/i 



uj - uj th + i«r - i/2r a 



(6.1) 



where <r„ is the complex coupling constant of the reso- 
nance to the n-th mode in the waveguide, T is the decay 
rate, and 1 /r a the amplification rate. The laser threshold 
is at Tr a = 1. 

We assume that the incide nt ra diation h as fr equ ency 
luq = uj th . Substitution of Eq. into Eq. (|4~T|) or ( flfl ) 
gives the simple result 



T 



-2/E 



2 A 



(6.2) 



for the limiting value of the noise figure on approaching 
the laser threshold. The limit is the same for detection 
in transmission and in reflection. Since the coupling con- 
tant |cr mo | 2 to the mode uiq of the incident radiation can 
be much sm aller than the total coupling constant E, the 
noise figure (6J!) has large fluctuations. We need to con- 
sider the statistical distribution p(J-) in the ensemble of 
random media. The typical (or modal) value of T is the 
value ^typ at which p(T) is maximal. We will see that 
this remains finite although the ensemble average -< T >- 
of T diverges. 



A. Waveguide geometry 

We first consider the case of an amplifying disordered 
waveguide. The total coupling constant E = E; + E r 
is the sum of the coupling constant E; = Yl n =i I "! 2 *° 
the left end of the waveguide and the coupling constant 



E 



2N 

n=N+l 



\cr n \ to the right. The assumption of 



equivalent channels implies that 



4/iV 



(6.3) 



Since the average of XfT is finite, it is reasonable to 



assume that 



typ 



i< 1/T -4/AT, or T x 



typ 



47V 



for complete population inversion. The scaling with N 
explains why the large- N theory of the previous section 
found a divergent noise figure at the laser threshold. We 
conclude that the divergency of T at L/£ a = ir in Fig. ^ 
is cut off at a value of order N, if T is identified with the 
typical value .Ftyp- 



B. Cavity geometry 

In the case of an amplifying disordered cavity, we can 
make a more precise statement on p(JF). Since there is 
only reflection there is only one E = E^Li \ a n\ 2 ■ The 
assumption of equivalent channels now gives 
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FIG. 4. Probability distribution of the noise figure near the 
laser thre shol d for an amplifying disordered cavity, computed 
from Eq. j6.7| ) for / = —1. The most probable value is T = N, 
while the average value diverges. 



2fN 



(6.4) 



Following the same reasoning as in the case of the wave- 
guide, we would conclude that .Ftyp ~~< I/-? 7 ^~ 1 = 
—2fN. We will see that this is correct within a factor of 
two. 

To compute p(JF) we need the distribution of the di- 
mensionless coupling constants u n — a n j\fT,. The N 
complex numbers u n form a vector u of length 1. Accord- 
ing to random-matrix theory |t]], the distribution p(S) of 
the scattering matrix is invariant under unitary transfor- 
mations S -> USU T (with U an N x N unitary matrix). 
It follows that the distribution p(u) of the vector u is 
invariant under rotations u — > Uu, hence 



p(u 1 ,U 2 , ■ ■ -,U N ) OC 6 ^1 - ^ \ u n\ 2 ^j 



(6.5) 



In other words, the vector u has the same distribution as 
a column of a matrix that is uniformly distributed in the 
unitary group ]Tq| . By integrating out N — 1 of the u n 's 
we find the marginal distribution of u mo , 



P(«m ) 



N- 1 



( 1 - Wr, 



,N-2 



(6.6) 



for TV > 2 and |u mo | 2 < 1. 

The distribution of T = — 2/|u mo | -2 becomes 

p(F) = -2f(N-l)(l + ?Pj N 2 T-\ 



(6.7) 



for N > 2 and T > ~2f. We have plotted p(F) in Fig. @ 
for complete population inversion (/ = — 1) and several 
choices of N. It is a broad distribution, all its moments 
are divergent. The typical value of the noise figure is the 
value at which p{T) becomes maximal, hence 
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ft 



typ 



-fN, N>2 



(6.8) 



In the single- mode case, in contrast, T = — 2/ for every 
member of the ensemble [hence p{T) — 8{T + 2/)]. We 
conclude that the typical value of the noise figure near 
the laser threshold of a disordered cavity is larger than 
in the single-mode case by a factor N/2. 



VII. ABSORBING MEDIA 

The general theory of Sec. || can also be applied to an 
abso rbin g medium, in equilibrium at temperature T > 0. 
Eq. (2.4) then has to be replaced with 



a out {u) = S(uj)a ia (u) + Q{u)b{uj) , (7.1) 
where the bosonic operator b has the expectation value 
(tf n (u)b m (u')) = 5 nm S(uj - u')f(cj,T) , (7.2) 
and the matrix Q is related to S by 

QQt = t-SSK (7.3) 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

FIG. 5. Excess noiS6 pOWGf Pex.c for an absorbing (solid 
line, left axis) respectively amplifying disordered waveguide 
(dashed line, right axis), in units of a 2 l\f\Io/L. The top panel 
is for detection in transmission, the bottom panel for detection 
in reflection. 




FIG. 6. Excess noise power P eX c for an absorbing (solid 
line, left axis) respectively amplifying disordered cavity 
(dashed line, right axis), in units of a 2 |/|7o. 



The formulas for F(£) of Sec. [II remain unchanged. 

Ensemble averages for absorbing systems follow from 
the corresponding results for amplifying systems by sub- 
stitution r a — > — r a . (This follows from a general duality 
theorem [|l7| between absorbing and amplifying systems.) 
The results for an absorbing disordered waveguide with 
detection in transmission are 



- Aal s 
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3L sinh s 
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2a 2 l 
3L 
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sinh s 



s cotanh s — 1 
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(7.4) 



(7.5) 



where s = L/^ a with £ a the absorption length. Similarly, 
for detection in reflection one has 



I = alo 
2a 2 / 



1 



41 
3L l 



-s cotanh s 



(7.6) 
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3L 



flos 



2 cotanh t 



1 



s cotanh s — 1 
sinh 3 s 



sinh s 

s 



cotanh s 
sinh 2 s 



sinh s 



(7.7) 



These formulas follow from Eqs. (5.1)-( p.4j ) upon substi- 
tution of s — > is. 

For an absorbing disor dered cavity, we find [substitut- 
ing 7 -> - 7 in Eqs. Q-(H)] , 



I = aln 



1 



1 + 7 



2a 2 // 7 



7 2 + 7 + 1 
(1 + 7) 4 



(7.8) 
(7.9) 



with 7 the dimensionless absorption rate. 

Since typically / -C 1 in absorbing systems, the noise 
figure T is dominated by shot noise, T ~ Iq/I. Instead 
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of T we therefore plot the excess noise power P oxc in 
Figs. H and [| In contrast to the monotonic increase 
of P exc with l/r a in amplifying systems, the absorbing 
systems show a maximum in P exc for certain geometries. 
The maximum occurs near L/£ a = 2 for the disordered 
waveguide with detection in transmission, and near 7 = 1 
for the disordered cavity. For larger absorption rates the 
excess noise power decreases because I becomes too small 
for appreciable beating with the spontaneous emission. 



VIII. CONCLUSION 

In summary, we have studied the photodetection 
statistics of coherent radiation that has been transmit- 
ted or reflected by an amplifying or absorbing random 
medium. The cumulant generating function P(£) is the 
sum of two terms. The first term is the contribution 
from spontaneous emission obtained in Ref. f§. The sec- 
ond term P exc is the excess noise due to beating of the 
coherent radiation with the spontaneous emission. Equa- 
tion ( |3.6| ) relates P exc to the transmission and reflection 
matrices of the medium. 

In the applications of our general result for the cu- 
mulant generating function we have concentrated on the 
second cumulant, which gives the spectral density P exc 
of the excess noise. We have found that P cxc increases 
monotonically with increasing amplification rate, while 
it has a maximum as a function of absorption rate in 
certain geometries. 

In amplifying systems we studied how the noise figure 
T increases on approaching the laser threshold. Near the 
laser threshold the noise figure shows large sample-to- 
sample fluctuations, such that its statistical distribution 
in an ensemble of random media has divergent first and 
higher moments. The most probable value of F is of 
the order of the number N of propagating modes in the 
medium, independent of material parameters such as the 
mean free path. It would be of interest to observe this 
universal limit in random lasers. 
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APPENDIX: DERIVATION OF EQ. (3.6) 



To ev aluate the Gaussian averages that lead to 



Eq. (3.6), it is convenient to use a matrix notation. We 
replace the summation in Eq. (3.2) by a multiplication 
of the vector a° ut with the projection Pa out , where the 



projection matrix V has zero elements except V n 
N + 1 <n< 2N. We thus write 

W = a [ dta out ^{t)Va out (t) . 
Jo 



(Al) 



Insertion of Eqs. (2.4) and (3.3) gives 



o Jo 



OO pQO 



W = j- \ dt I du / du/[a int MS t M + cHv^M] 



xP [S(w> ir V) + V(lj')c\lj')] 



(A2) 



As explained in Sec. Ill we discretise the frequency as 
Lu p = pA, p = 1,2, 3,.... The integral over frequency is 
then replaced with a summation, 



Jo 



P =i 



We write Eq. (A.2) as a matrix multiplication, 
CW = a in Ua in + cPct + a^CV + cCa in 
with the definitions 



(A3) 



(A4) 



.4 



np,n'p' 



2tt 



2^ Jo 



B np , nlp/ =^ I dt {v\u> v )VV{u p ,)). e^b-rt* 



= A 1 /^^) • (A5) 



We now apply the optical equivalence theorem [ fTl) , 
as discussed in Sec. HI. The operators a™ p are replaced 
by constant numbers S nTno S PPo (2TTlo/ A) 1 / 2 . The opera- 
tors c np are replaced by independent Gaussian variables, 
such that the expectation value ([O]) takes the form of a 
Gaussian integral, 



5 iW :) = J rf{c„ p }exp 



np 

= J d{c np } exp [a in *Aa in - cMc* + a in *C^c* + cCa in ] , 

(A6) 

where we have defined 



$nn' $pp' 



(A7) 



We eliminate the cross-terms of a m and c in Eq. ( |Aq ) 
by the substitution 



c' = c* - M~ 1 Ca in , 



(A8) 



leading to 



7 



(: e 6w :) = exp [a in *(A + C t Af- 1 C7)o ln ] 

x J d{c' np }exp(-c'Mc'*) . (A9) 

The integral is proportional to the determinant of M _1 , 
giving the generating function 

= constant - In ||M|| + a in * (A + CftM^C) a in 
= constant — In ||M|| 



+ ^( A + cHi- 1 c) 



mapa.mapa 



(A10) 



The additive constant follows from F(0) = 0. The term 
— In || M || is the contribution from amplified spontaneous 
emission calculated in Ref. || . The term proportional to 
Io is the excess noise of the coherent radiation, termed 
F rxr in Sec. [II. 



Eq. (A1C) can be simplified in the long-time regime, 
uj c t ^S> 1. We may then set A = 2-k/t and use 



e i&(p- P ')t dt _ r6pp/ 



(All) 



The matrices defined in Eq. (A5) thus become diagonal 
in the frequency index, 



.4 



np,n'p f 



2tt 



(S\uj p )TS(co p )) nn ,6 ppl (A12) 



and similarly for B and C. We then find 

(A + C Ad 1 C) np ^ n i p i — 

n ^ T (&V[l + atfWiV]- 1 S) nn ,6 ) 



2tt 



<O pp ' , 



(A13) 



where /, S, an d V are evaluated at m — lu p . Substitution 
into Eq. flAlOD gives the result tfHfy for F oxc (£). 

Simplification of Eq. ( |Alfj| ) is also possible in the short- 
time regime, when J7 c r <C 1, with f2 c the frequency range 
over which SS^ differs appreciably from the unit matrix. 
The generating function then is 



~2T 



dwf(w,T) 



x (1 -t{uj)t\uj) -t{uj)t\uj)) t(u )) ■ (A14) 
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